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Abstract 

As we all known, there is still a long way for us to solve arbitrary mul- 
tivariate Lagrange interpolation in theory. Nevertheless, it is well accepted 
that theories about Lagrange interpolation on special point sets should cast 
important lights on the general solution. In this paper, we propose a new 
type of bivariate point sets, quasi-tower sets, whose geometry is more natural 
than some known point sets such as cartesian sets and tower sets. For bivari- 
ate Lagrange interpolation on quasi-tower sets, we construct the associated 
degree reducing interpolation monomial and Newton bases w.r.t. common 
monomial orderings theoretically. Moreover, by inputting these bases into 
Buchberger-Moller algorithm, we obtain the reduced Grobner bases for van- 
ishing ideals of quasi-tower sets much more efficiently than before. 

Key words: Bivariate Lagrange interpolation. Quasi- Tower set. Degree 
reducing interpolation basis, Buchberger-Moller algorithm 
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1. Introduction 

Given a set H = {^^^\ . . .,^^''^} C of /i distinct points with F a field. 
For prescribed values G F, i = 1, . . . , /i, bivariate Lagrange interpolation is 
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to find all polynomials p satisfying 

p{^^'^) = fi, ^ = l,...,/i, (1) 

where 11^ := F[x,?/] is the bivariate polynomial ring over F. 

In most cases, for fixed monomial ordering -<, only the bases for the 
degree reducing interpolation space V C 11^ (see [jj) are pursued for Lagrange 
interpolation ([1]). Let {pi, . . . be a basis for V satisfying pi ^ p2 -< ■ ■ ■ -< 

Pm- If 

for certain reordering of S, then we call {pi, . . . a degree reducing in- 
terpolation Newton feasis(DRINB) w.r.t. -< for ([1]). Moreover, A/'-^(H), the 
Grobner escalier of X(H) w.r.t. -<(see jif), is the degree reducing interpolation 
monomial 6as^s(DRIMB) w.r.t. -< for ([1]). 

Let X(H) C be the vanishing ideal of H and (H) the reduced Grobner 
basis for X(H) w.r.t. -<. In 1982, jsf proposed Buchberger-Moller(BM for 
short) algorithm for computing bases of vanishing ideals. Specifically, with 
inputted S and -<, BM algorithm outputs ^^(H), A/'^(S), and a DRINB for 
Lagrange interpolation ([1]), namely arbitrary bivariate Lagrange interpola- 
tion can be solved with it algorithmly. Nonetheless, the poor complexity of 
BM algorithm(see severely constrains its applications. 

Since the geometry of H is a dominant factor for the structures of (5) 
and A/'^(H), theories about Lagrange interpolation on point sets with describ- 
able special geometries will certainly throw light on the general theoretical 
solution of ([1]). jsilGf studied multivariate Lagrange interpolation on cartesian 
sets(aka lower point sets) and investigated bivariate Lagrange interpola- 
tion on tower sets. 

This paper will introduce, in Section [3l a new type of bivariate point sets, 
quasi-tower sets, that is more natural than cartesian and tower sets. As the 
section of main results. Section |4] will discover the degree reducing interpo- 
lation bases w.r.t. certain common monomial ordering such as graded lex 
order for bivariate Lagrange interpolation on quasi-tower sets theoretically. 
Finally, Section[5]will give an algorithm for computing reduced Grobner bases 
for vanishing ideals of quasi-tower sets. 

The following section will serve as a preparation for the paper. 
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2. Preliminary 

Let No and be the monoid of nonnegative integers and bivariate mono- 
mials in respectively. A monomial in has the form X" := x^^y"'^ for 

CX = («!, Oo) G Nq. 

As in |8[, we use -<icx, -<inviex5 -<griex5 and -<grcvicx to represent lexicographic 
order, inverse lexicographic order, graded lex order, and graded reverse lex 
order respectively. Hereafter, we always assume that y -< x for any fixed 
monomial ordering -<. 

For a nonzero / G H^, we let LM^(/) and S^{f) signify the leading 
monomial and leading bidegree(LM(f) = X^^^^^^) of / w.r.t. -< respectively. 
Furthermore, set 

f^g:= 6M)^S^i9)- 
Let A be an arbitrary finite subset of Nq. We say that it is lower if 

R(a) := {(«;, a'2) e Nl : < a'i < ai,i = 1,2} C A 

always holds for any a = (01,02) G A. When A is lower, we let nij = 
max(/i ,,)g_4 h,0 < j < u, with u = max(o,fc)g^ k. Since A can be determined 
uniquely by {u + l)-tuple (mo, mi, . . . , rriu), we represent it as 

A = L^{mo,mi, . . . ,my). 

Likewise, it can also be represented as Ly(no, . . . , rimo) with ni = max(j^fc)g_4 k, 
< i < rrin. 

As in j7|, we constructed two particular lower sets SxC^), Sy{E) C Ng 
from S. In more details, we cover H by lines /g, , . . . , parallel to the x- 
axis and assume that, without loss of generality, there are rrij + 1 points, say 
Uqj, u^j, . . . , u^^.J, on Ij with mo > mi > ■ ■ ■ > m,^ > 0. Set 

SxiE) := : < z < m„ < J < u}, 

which equals to L2:(mo, mi, . . . , m,y). In like fashion, we cover S by lines 
/q, /i, . . . , parallel to the y-axis and denote the points on Zf by u^q, . . . , „. 
with no > ■ ■ ■ > nx > 0- We have 

Sy(E) := : < z < A, < j < ni} = Ly{no,ni, . . .,nx). 

In addition, the sets of abscissae and ordinates are defined as 

H,{E) ■.={x : {x, y) G /J n H}, < j < v, 
V^,(S) :={y : G/fnS}, < z < A. 
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A point set S satisfying S^C^) = SyiS) is called a cartesian set that has 
the following property: 

Proposition 1. ^ point set S C is cartesian if and only if 

Ho{E) D H,{E) D ■ ■ ■ D H,{E) 

or 

VoiE)DV,{E)D---DVx{E). 

Enlightened by the notion of cartesian sets, 0] introduced tower sets in 
as follows. 

Definition 1. j3] We say that a set S of distinct points in is x-tower if 
lower set Sx{E) = Lj:{mo, mi, . . . , m,^) C Ng with mo > mi > ■ ■ ■ > m,^ > 
such that 

E := {{xij,yj) : (z,j) G S^iE)}, 

where Xij G Hq{E), {i,j) G Sx{E), are distinct for fixed j. Its full name is 
Sx(E)-x-tower set. Similarly, if lower set Sy{E) = Ly{no,ni, . . . ,nx) C Ng, 
no > ni > ■ ■ ■ > nx > 0, such that 

E := {{x^,y^j) : G 53,(H)}, 

where yij G V()(S),(i,j) G Sy{E), are distinct for fixed i, we will call H a 
y-tower set (or Sy{E)-y-tower set in its full name). 

About the bivariate Lagrange interpolation on a tower set, Q proved the 
succeeding theorem. 

Theorem 2. (3] Given an x-tower set H C F^. The DRIMB w.r.t. -<griex or 
-<icx for ^ is 

{xY ■■{^,3)^S,{E)}. 
IfE is y-tower, then the DRIMB w.r. t. ~^grcvlGx or ~>invlcx 

for is 

{xy ■.{z,j)eSy{E)}. 
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3. Quasi-Tower Sets 



In this section, we will introduce the notion of generalized tower and 
quasi-tower sets and compare them with tower sets and cartesian sets. 

Definition 2. Let S be a set of distinct points in F^. We call S a generalized 
x-tower set if S^C^) = Lx{mQ,mi, . . . ,m^) C Nq with mg > mi > ■ ■ ■ > 
TTiiy > such that 



where distinct for fixed j. In like fashion, if lower set 5*^(2) = 

Ly(no) ni, . . . , n\) C Ng, tiq > rii > • ■ ■ > n\ > Q , such that 



where yij are distinct for fixed i, we will call S a generalized y-tower set. 

Observe Definition [T] and [21 We find that they are similar to each other 
except two more conditions in Definition [T] which implies that a tower set 
is always generalized tower, namely the notion of generalized tower sets is 
actually a generalization of tower sets'. Therefore, what interests us most 
now should be the generalized tower sets that are not tower. 

Definition 3. If H C is a generalized x-tower set that is not x-tower, then 
we call it a quasi-x -tower set. Similarly, a non-y-tower generalized y-tower 
set is called a quasi-y -tower set. 

From Definition |2] and |3] we know that the geometry of quasi-tower sets 
is much freer than tower sets' in the sense that they need not fulfill the two 
extra conditions that demand Xij C HqCE) for x-tower and Xij C Vo(H) for 
l/-tower. The following propositions will show us how did quasi-tower sets 
acquire their name. 

Proposition 3. Let S be a quasi-x-tower set with S^C^) = Lrc{mQ, . . . ^m^). 
If B G is a set of distinct horizontal points such that i? fl H = and 
[J'^^^Q Hi{E) C Hq{B), then EU B is x-tower and will be known as a derived 
x-tower set from E with base set B. 



{i,j)eSy{E)}, 
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Proof. By hypothesis, we assume that Sy(E) = Ly{nQ,ni, . . . ,nx) and 
B = {(x^,?/^), (xf ,?/^), . . . , (x^^,?/^)} where xf ^ xf,i 7^ j, and ^ 
[j^=oVj{E). Recall the construction process of S's(S) in Section |2l Since 
[jj^^Q Hi{E) C Ho{B) and S is quasi- x-tower, we can deduce easily that 

Sx(E U B) = L^{mB, mo, ... , m^), 

where > mg > ■ ■ ■ > m,^ > 0. From Definition [2] and [31 S is quasi- x-tower 
implies that 

where x^- are distinct for fixed j, hence, we have 

SU5 = {(xij,%) : G S^{EUB)} 

where 

Xio = xf,yo = y^, i = 0, 1, . . . , m^, 

Xjj+i = xj^., i/j+i = yj, 2 = 0, 1, ... , mj, j = 0, 1, . . . , z/, 

and Xjj are distinct for fixed j. Therefore, by Definition [H S U -B is x- 
tower. □ 

In the same way, we can prove the following proposition for quasi-y-tower 
sets. 

Proposition 4. Let E be a quasi-y-tower set with 5*^(2) = L,y(no, rii, . . . , nx) ■ 
If B G is a set of distinct vertical points such that i? fl H = and 
U;=o^.(S) C Vo{B), thenEUB IS y-tower and will be named a derived 
y-tower set from E with base set B. 

Proposition [3] and H] imply the following corollary immediately. 

Corollary 5. Let E g¥^ be a quasi-tower set. Then there are infinite derived 
tower sets from E. 

Example 1. 

E = {(0.2, 0.4), (0.4, 0.4), (0.8, 0.4), (1, 0.4), (1.2, 0.4), (1.6, 0.4), 
(2, 0.4), (0, 0.6), (0.6, 0.6), (1.4, 0.6), (1.8, 0.6), (2.2, 0.6), 
(0.6, 1.1), (1, 1.1), (1.2, 1.1), (1.8, 1.1), (1.2, 1.45), (1.8, 1.45), 
(2.4, 1.45), (0.2, 1.8), (1.4, 1.8)} C Q2 
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is a quasi- x-tower set (shown in (a) of Figured]) and 

5 = {(0, 0), (0.2, 0), (0.4, 0), (0.6, 0), (0.8, 0), (1, 0), 

(1.2, 0), (1.4, 0), (1.6, 0), (1.8, 0), (2, 0), (2.2, 0), (2.4, 0)}. 

Obviously, H U i?(shown in (b) of Figure [1]) is an x-tower set that is derived 
from S with base set B. 
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(a) S (b) S U B 

Figure 1: A quasi- x-tower set S and a derived a;-tower set SUB from S with base set B 

Proposition 6 of 0] discovered the relation between cartesian sets and 
tower sets. The next proposition shows the relation between cartesian sets 
and quasi-tower sets. 

Proposition 6. A quasi-tower set is not a cartesian set and vice versa. 

Proof. From Definition [3] we can deduce easily that a point set S C is 
quasi- x-tower if and only if it is a generalized x-tower set such that there 
exists at least one point ^ G S whose abscissa is not in HqCE). Proposition [T] 
implies the statement of this proposition immediately for quasi- x-tower sets. 
The quasi-y-tower cases can be proved similarly. □ 

4. Bases for Lagrange Interpolation on Quasi-Tower Sets 

In this section, we will present the degree reducing interpolation bases for 
Lagrange interpolation on quasi-tower sets w.r.t. -<icx, -<grcvicx etc. theoreti- 
cally. Let us begin with the following lemmas. 
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Lemma 7. Fix monomial ordering -<gricx- Let 

E' = {{Xij,yj) : (ij) G S^{E')} 

be an x-tower set in with SxiE!) = L^.(mo,mi, . . . ,my). If = 
{91, 92,---, gs} with Qs -<gricx " " " -<gricx 91, then 

(i) s = u + 2. 

(ii) gi = {x - xofi){x - xifi) ■ ■ ■ (x - x^o.o)- 

(iii) 9s = {y- yo){y -yi)---{y- yu)- 

(iv) 9j e (y - yo) and 

LM^,J9,) = ^M^,.,Jg,) = x^^-^+'y^-\ j = 2, . . . , s - 1. 

(v) X({(xo,o, yo), • • • , (a^mo.o, yo)}) = {9i,y- Vo)- 

Proof. From Theorem|2]we know that A/'-^g^j^^(S') = {x'y^ : (z,j) G S^iE')}. 
Since Sx{E') = L^(mo, mi, ... , m^), it is easily checked that 

{LM^^,^^(^?,) ■.l<z<s} = {x'^^+'y^ ■.0<j<u}U {y^^'}, 

which imphes (i) readily. In the proof of Theorem |2]^For brevity, we refer the 
reader to (3] for details), we find that {x — Xofi){x — xi^) ■ ■ ■ {x — Xmofl), {y — 
yo){y - yi) ■ ■ ■ {y - yu) e Q^^^^^^{E'). Since S' is x-tower, we have mo > mi > 
■ ■ ■ > ruv > 0. Thus, for < j < z/ — 1, 

degx^'^^^y^ = mj + l+ j 

> mj+i + 1 + j + 1 
= degx"'^+'+^y^+\ 

Next, we claim that deg x^^J^^y^ > deg y"'^^, that is m^ > z/ — j, j = 0, . . . , u. 
In fact, it is easy to see from mo > mi > ■ ■ ■ > m,^ > that ruj — my > v —j, 
hence, m^ > z/ — j, j = 0, . . . , i^. Therefore, (ii), (iii) follows. 

Still from the proof of Theorem |2l we know that gj G — yo) and there 
is no monomial x'^^y"'' of gj except LM^^^^^^{gj) satisfying ai > rnj-i + 1. 
Therefore, the definition of -<iex implies (iv). 

Finally, (v) is a direct consequence of Lemma 8 in [3] . □ 

Lemma 8. (sj Let E and E' be varieties in W^. Then 

I{E):I{E')=I{E-E'). 
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Lemma 9. [sl Let X he an ideal and g an element of IP . If {hi, . . . , hp} is 
a basis of the ideal X fl {g) , then {hi/g, . . . , hp/g} is a basis ofX : (g). 

The following theorem shows us the degree reducing interpolation bases 
for bivariate Lagrange interpolation on quasi- x-tower sets w.r.t. -<griex- 

Theorem 10. Let S be a quasi-x-tower set in F^. The DRIMB for Lagrange 
interpolation on S w.r.t. -<gricx is 

AA^^,^^(S) = {xy:(z,j)e5.(S)}. 

Proof. Firstly, we assume that S^i^) = L^(mi, . . . , m;^). By Corollary 
[5l there exists a derived x-tower set from H, say S', with base set B = 
{(zo, yo), (xi, yo), • • • , {xmo, Uo)}- Thus, Proposition [3] implies that S^{E') = 
La;(^05'"^i; • • • I'^v)- Moreover, it follows from (v) of Lemma [7| that 

X{B) = {{x - xo){x - xi) ■ ■ ■ (x - XmJ, y - yo) 

= ((X - Xo)(x - Xi) ■ ■ ■ (X - Xmo)) + {y- yo)- 

Since 5 n S = 0, Lemma E leads to X(S) = X(S' - B) = J(S') : X{B). 
Therefore, we have 

X(H) = X{E') : (((x - xo)(x - xi) ■ ■ ■ (x - x„J) + {y - yo)) 

= (X(H') : ((x-xo)(x-xi)---(x-x„J)) H (X(H') : (y-yo)). 

To be sure, it follows from (ii) of Lemma [7] that ((x — xo)(x — xi) • ■ ■ (x — 
Xmo)) C X(S'), hence, 

X(H') : ((x - Xo)(x - xi) ■ ■ ■ (x - Xmo)) = 

which means 

X(H)=X(H'):(y-yo). 

Recalling Lemma [TJ we suppose that = {gi, g2, ■ ■ ■ , gu+2} with 

gu+2 -<gricx ■ ■ ■ -<gricx gi- Hcuce, wc havc 

gi = {X - Xo)(x - Xi) ■ ■ ■ (X - Xmo), 

9i G {y-yo), i = 2,...,u + 2. 
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Next, for the basis for X(S), we should compute the basis for X(S') n{y — Uq) 
so that we can apply Lemma O Let t be a new variable. According to [sl, 

X(S') n (y - y,) = (tJ(S'), (1 -t){y- yo)) D 

= {tgi,tg2,...,tg^+2,i^ -t)iy-yo)) ntf. 

Since gt e {y - yo), i = 2, . . . , u + 2, gi/{y - yo) E H^, hence, 

tgi + {gi/{y - 2/0)) (i -t){y- yo) 
={gi/{y -yo))t{y ~yo) + {gi/{y - yo)){^ - t){y - yo) 
={gi/{y -yo)){y - yo) 
=gi, 

holds for i = 2, . . . ,1/ + 2. Thus, 

X(H') n {y - yo) = {tg,, (1 - t){y - yo),g2, ■ ■ ■ , gu+2) n U\ 

Since Q^^^^^^{E') = {gi,g2, ■ ■ • ,fi'v+2}, it is easily checked from Buchberger's 
S-pair criterion and (iv) of Lemma [7] that 

{tgi, (1 - t){y - yo),g2, ■ ■ ■ , gu+2} 

is the reduced Grobner basis for (tX(S'), (1 — t){y — yo)) w.r.t. -<iex with 
y -<icx X -<icx t. Consequently, 

X(S') n{y- yo) = (^2, • • • , gu+2), 

hence, 

I{E) =X(S') : (y-yo) = (^2/(2/ - 2/o), ^^+2/(2/ - 2/o)), 
follows from Lemma O In fact, we can deduce that 

^^g.icx(^) = {g2/{y-yo),---,gu+2/{y-yo)} 
=■■ {g'i,---,g'u+i}- 

By Lemma [71 S' is x-tower implies that 

{LM^^^,^^(^7,) : 1 < J < z/ + 2} = {x'^^-^'y^ : < j < z/} U {y-'^'}, 
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therefore, 

{LM^^,^^(^7;) : 1 < J < z/ + 1} = : I < J < i^} U {y^}. 

As a result, since S^i^) = L^{mi, . . . , m^), we have 

which finishes the proof. □ 

Theorem [TO] has solved the DRIMB problems w.r.t. -<griox- Next, let us 
turn to DRINB problems. 

Theorem 11. Let S C be a quasi-x-tower set of points 

which give rise to polynomials 

j-i i-l 



^"j Yliv - yt) - x,j), (i, j) e ^,(2), 



t=0 s=0 



where (p^^ = ^/Y[t=oiyj ~ Vt) ns=o('^«i ~ -^sj) ^ ^' ^'^'^ empty products 
are taken as 1. Then 

is a DRINB w.r.t. -<gricx for satisfying 

^U'^mn) = %i),(m,n), ih j) ^invlex {m,n). 

Proof. Fix {i,j) G S^i^)- When = (m, n), since Xqj ^ Xij ^ ■ ■ • ^ Xij 
and I/O 7^ Z/i 7^ ■ ■ ■ 7^ 



i-i j-i 



t=0 s=0 



follows. Otherwise, if (i, j) ^mvicx ijn^n), we have j > n, or j = n,i > m. 
When j > n, we have 



s=0 
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and when j = n,i > m, 
i-i 

t=0 
n-1 

'/'jj J_ J_ (.Vn Vt) i,'^mn •^On) ' ' ' ij^mn •^mn) ' ' ' ij^mn ■^i—l,ri) 
t=0 

= 0, 

which leads to 

0fj(Wmn) = 0, ^invlcx (m,n), 

namely Qx is a Newton basis for Span^Qz- By Theorem [T0|, it is easy to see 
that SpanpQz = SpanpA/".^ riex('^)- Therefore, is a DRINB w.r.t. -<griex for 

m- □ 

Similarly, we can prove the following theorems: 

Theorem 12. Let S be a quasi-y -tower set in F^. The DRIMB for Lagrange 
interpolation on S w.r.t. -<grcvicx is 



Theorem 13. Let Ec¥^ be a quasi-y -tower set of points 

We define the polynomials 

i-i j-i 

s=0 t=0 

where ip^j = l/ns=o(^« ~ •^«) nt=o(y«j ~ Hit) ^ The empty products are 
taken as 1. Then, 

Qy = {<|>'^,■.{^,3)eSyiE)} (2) 
is a DRINB w.r.t. -<grcvicx for ([I]) satisfying 
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So far, we have discovered the degree reducing interpolation bases for 
Lagrange interpolation on quasi-tower sets w.r.t. -<griGx and -<greviex- Next, 
we should turn to -<icx and -<invicx cases. 

Lemma 14. ^ Let 

= {«mn = i^mn^ymn) ■ {^,1%) G Sy{E)} 

be a set of distinct points in F^. Then 

(i) the set {xV : ihj) e is the DRIMB as well as {(pf^ : e 
Sx{E)} is a DRINB w.r.t. -<icx for (P, where 

i-i i-l 

^% = U^y - ylt) - {^.J) e ^.(H), 

t=0 s=0 

OTi/i (/pfj = l/nt=o(?/oj ~ ns=o(^ij' ~ ^sj) ^ -'^ '^'^'^ empty products 
taken as 1; 

(ii) the set {xV : e ^^(H)} 2S t/ie DRIMB as well as {^^ : (ij) E 
5*^(2)} is a DRINB w.r.t. -<invicx for ([I]), where 

^% = - ^'o) Uiy - yl), ihj) 6 SyiE), 

s=0 t=0 

with (f\j = 1/ ns=o(^io ~^so) ntol^ij- ~yit) ^ ^^'^ empty products are 
taken as 1. 

Finally, Theorem [T01IT3] and Lemma [E] together give rise to our main 
theorem. 

Theorem 15. Let 

he a quasi-x-tower set in F^. Then the set {x^y^ : G Sx(E^)} is the 

DRIMB as well as {(p^j : {i,j) G Sx{E^)} is a DRINB w.r.t. ^lex or ^griex 
for ([1]), where 

i=0 s=0 
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with If = ^/ Y[i=oiyj ~yt) YYsJoi^ij~^sj) ^ ^ ^'^'^ empty products taken 
as 1. 

If 

^' = {ul = {x\,yl)■.{^,3)eSy{Ey)}(l¥' 

is a quasi-y -tower set. Then the set {x^y^ : {i,j) G S'y(S^)} is the DRIMB as 
well as {(p^ij : G 5*^(2^)} is a DRINB w.r.t. -<invicx or -<gi.cvicx for (P, 

where 

wt/i v^f^- = 1/ ns=o(^f ^s) ntol^ij- ~ ^il) ^ "'^^ empty products are 
taken as 1. 

Example 2. Let us continue with Example [H Observeing Figure [H we have 

=1,(6,4,3,2,1). 
Therefore, according to Theorem [T5l 

={1, X, x^, a;^ a;^, a;^ a;^ y, xy, x^y, x^y, x^y, (3) 

2 222 S 24 424^ 

y ,xy ,x y ,x y ,y ,xy ,x y ,y ,xy } 
and a DRINB for interpolation on S w.r.t. -<icx or -<griex is 

25 

1, 5(x - 0.2), —(a; - 0.2)(x - 0.4), . . . , 
6 

25 

5(y-0.4),yx(y-0.4),..., (4) 

20, , 1 

y(2/-0.4)(y-0.6),... L 

5. Algorithms and Timings 

At the beginning of this section is a restatement of the classical BM 
algorithm with the notation established above. 
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Algorithm 1. (BM Algorithm) 

Input: A set of distinct points S = {^^^^ : i = 1, . . . , /i} C F'^ and a fixed 
monomial ordering -<. 

Output: The 3-tuple {G,N,Q), where G is the reduced Grobner basis 
for X(H) w.r.t. -<, N is the Grobner escalier of X(H) (the DRIMB for 1^ 
also) w.r.t. -<, and Q is a DRINB w.r.t. -< for dU). 

BMl. Start with lists G = [],N = [],Q = [],L = [1], and a matrix 
B = (bij) over F with fi columns and zero rows initially. 

BM2. If L = [ ], return {G,N,Q) and stop. Otherwise, choose the 
monomial t = min^L, and delete t from L. 

BM3. Compute the evaluation vector {t{^^^^), . . . ,t{^^^^)), and reduce it 
against the rows of B to obtain 

i 

BM4.. If (t>i, . . . , v^) = (0, . . . , 0), then append the polynomial t—"^. aiQi 
to the list G, where qi is the ith element of Q. Remove from L all the multiples 
of t. Continue with BM2. 

BM5. Otherwise {vi, . . . , Vfj) ^ (0, . . . , 0), add {v\, . . . , v^) as a new row 
to B and t — a^qi as a new element to Q. Append the monomial t to A, 
and add to L those elements of {x\t^ . . . , x^t} that are neither multiples of 
an element of L nor of LM^(G) := \\M^{g) : g^G}. Continue with BM2. 

Like TBM algorithm in 0], we have the following QTBM algorithm for 
quasi-tower sets. 

Algorithm 2. (QTBM algorithm) 

Input: A quasi- a;-tower(quasi-?/-tower) set S C F^ of /i points and a fixed 
monomial ordering -<griex(-<greviex) or -<icx(-<mviex)- 

Output: The 3-tuple {G,N,Q), where G is the reduced Grobner basis 
for X(S), N is the Grobner escalier of X(S), and Q is a DRINB for ([1]). 

QTBMl. Construct lower set S,^{E){Sy{E)). 

QTBM2. Compute the sets N and Q according to Theorem [TH 

QTBM3. Compute the set L := {x ■ t : t e N} [j{y ■ t : t G A} \ A. 

QTBM4. Construct fix fi matrix G whose {h, k) entry is <Ph{'^k){4>h{'^\)) 
where 0^(0^), ul{ul) are hth and kth elements of Q = {4>ij{(p^j) '■ G 
S,{E){Sy{E))} and H = {<J<J : {m,n) E S.{E){Sy{E))} w.r.t. the 
increasing -<inviex("<iex) on and {m,n) respectively. 

QTBM5. Goto BM2 with A, Q, L, G for the reduced Grobner basis G. 
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Example 3. We continue with Example [T] and [21 Example [2] had obtained 
and Q in and (HD respectively. According to QTBM3 and QTBM4, 

we get 

L = {x^ x^y, x^y, x^y^ ^/^ xy^} 



and 



C 



( 1 









1 

3 
1 



Finally, QTBM5 yields 

G = {(y-0.4)(y 



0.6)(y-l.l)(2/-1.45)(y-1.8),...}. 



At the end, we will show the timings for the computations of BM-problems 
on quasi-tower sets in finite prime fields with size q w.r.t. monomial 
orderings -<griex and -<inviex respectively. QTBM and BM algorithms were 
implemented on Maple 12 installed on a laptop with 768 Mb RAM and 1.5 
GHz CPU. 

For field F37 and -<gricx5 



Algorithms 


#2 


300 


500 


800 


1200 




QTBM 




2.563 s 


8.773 s 


26.068 s 


63.772 s 


BM 




20.690 s 


80.375 s 


288.364 s 


901.246 s 



For field F43 and -<inviex, 



Algorithms 


#2 


500 


800 


1000 


1200 




QTBM 




8.171 s 


25.607 s 


48.170 s 


70.802 s 


BM 




114.996 s 


375.610 s 


881.988 s 


1347.287 
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